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The article is a natural continuation of our paper Quantum scalar field 

in FRW Universe with constant electromagnetic background, Int. J. Mod. 
Phys. A12, 4837 (1997). We generalize the latter consideration to the case 
of massive spinor field, which is placed in FRW Universe of special type with 
a constant electromagnetic field. To this end special sets of exact solutions 
of Dirac equation in the background under consideration are constructed and 
classified. Using these solutions representations for out-in, in-in, and out-out 
spinor Green functions are explicitly constructed as proper-time integrals over 
the corresponding contours in complex proper-time plane. The vacuum-to- 
vacuum transition amplitude and number of created particles are found and 
vacuum instability is discussed. The mean values of the current and energy- 
momentum tensor are evaluated, and different approximations for them are 
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presented. The back reaction related to particle creation and to the polariza- 
tion of the unstable vacuum is estimated in different regimes. 
PACS 04.62.+V, 12.20.Ds 

I. INTRODUCTION 

The article is a natural continuation of our paper Here we generalize the latter 
consideration to the case of a massive spinor field which is placed in Friedmann-Robertson- 
Walker (FRW) Universe of special type with constant electromagnetic field. First, we ought 
to repeat briefly a motivation for such an activity. 

It is quite well known fact that quantum field theory in an external background is, 
generally speaking, theory with unstable vacuum. The vacuum instability leads to many 
interesting features, among which particle creation from the vacuum is one of the most 
beautiful non-perturbative phenomenon. One has to treat it exactly in regard to the external 
field. The latter has been realized long ago by Schwinger 0. The particle creation effects 
together with back reaction issue are important in black hole physics and in dynamics of 
the early Universe (EU) (see, for example, and references therein). 

In quantum field theory with unstable vacuum it is necessary to construct different kinds 
of Green functions (GF), e.g. besides the causal GF (out-in GF) one has to use so called 
in-in GF, out-out GF, and so on [^j-^] (for a review and technical details see |^). General 
methods of such GF construction in electromagnetic background have been developed in 
IPJ^. A possible generalization of the formalism to an external gravitational background 
has been given in Ref. 0. Since Ref. it has been known that causal (out-in) GF may 
be presented as a proper-time integral over a real infinite contour. At the same time, in 
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the unstable vacuum case the in-in and out-out GF differ from the causal one. It was 
shown |IO| , p!T| there are examples of external fileds (electromagnetic with constant uniform 
invariants) when these functions may be presented by the same proper-time integrals (with 
the same integrand) but over another contours in the complex proper-time plane. Then, 
it is not difficult to compare contributions from the in-in GF and from the causal one. 
The complete set of GF mentioned is necessary for the construction of the Furry picture 
in interacting theories, and even in noninteracting cases one has to use them to define, for 
example, the back reaction of particles created and to construct different kinds of effective 



actions (EA) (for a general introduction to EA in background field method, see |T2|, and for 
review of modern generalizations, see ||13[). The such proper-time representation of GF may 
be the necessary step in the study of chiral symmetry breaking in QED and the four-fermion 
models under the action of gravitational and electromagnetic fields (see [0] and references 
therein). /^From another point, in-in GF which gives the origin to in-in EA maybe used in 
more realistic theories, like GUT theories with scalars, spinors and vectors in order to analyse 
the properties of above electro-gravitational background in the EU. For example, one of 
extremely interesting questions there is: can we realise the asymptotic conformal invariance 
phenomenon (which means that theory becomes approximately conformally invariant at 
large curvature) [|l^ even for in-in EA, or in other words for mean values in EU. Taking into 
account our recent study of in-in GF structure for scalars it looks quite interesting next 
application of above calculation. 

It may be likely that EU is filled with some type of electromagnetic fields. For example, 
recently (see ||T6|JT7[| and references therein) the possibility of existence and role of primordial 
magnetic fields in EU have been discussed. From another point the possibility of existence 
of electromagnetic field in the EU has been discussed long ago in |T8|JT9| . It has been shown 
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there that the presence of the electrical field in the EU increases significantly the gravitational 
particle creation from the vacuum. In principle, this process may be considered as a source 
for the dominant part of the Universe mass. 

Bearing in mind the above cosmological motivations it is becoming interesting to study 
the quantum field theory in curved background with electromagnetic field (of a special form 
in order to solve the problem analytically). In the present paper we are going to consider 
a massive spinor field placed in the expanding FRW Universe with the scale factor fl{r)) 
(in terms of the conformal time) ^^^(^7) = b'^rj^ + a^. Such a scale factor corresponds to the 
expanding radiation-dominated FRW Universe. In terms of physical time t the corresponding 
metric may be written as follows: 

ds^ = dt'^ -Q\t){dx'^ + dy^ + dz^) , (1) 

where for small times \t\ <^ a^/b, Q'^{t) ^ a^[l + (bt/a^)'^], and for large times \t\ S> 
a^/b, r2^(t) ~ 2b\t\ (see [|15]). Moreover, such FRW Universe will be filled by the constant 



electromagnetic field. 

Thus, we start from the theory of massive spinor in above background. Making a confor- 
mal transformation we remain with QED in fiat background but with time-dependent mass 
(QED-fi theory). In the Sect. II special sets of exact solutions of Dirac equation in QED- 
Q theory are constructed and classified as corresponding to particles and antiparticles at 
t — > ±00. In the Sect. Ill, using these solutions, representations for out-in, in-in, and out-out 
spinor GF are explicitly constructed as proper-time integrals over the corresponding con- 
tours in complex proper-time plane. As far as we know, it is a first explicit example for the 
proper-time representations for complete set of spinor GF in gravitational-electromagnetic 
background. In the Sect. IV we are interested in to reveal global features of the theory. The 
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vacuum-to-vacuum transition amplitudes and number of created particles are found and 
vacuum instability is discussed. It is seen the creation process is a coherent effect of both 
fields. The all mean values of the current and energy-momentum tensor are presented in 
the same manner as the proper-time integrals, and evaluated. The different approximations 
for them are investigated. The back reaction produced by both of particles created from a 
vacuum and polarization of an unstable vacuum estimated in different regimes. It is shown 
a behaviour of such components in time are quite different. 



II. CLASSIFIED SETS OF EXACT SOLUTIONS 

In this Section we study exact solutions of the Dirac equation in an external constant 
uniform electromagnetic background and in a time- dependent mass-like potential, which 
effectively reproduces effects of a gravitational background (solutions of the Dirac equation 
of QED-fi theory), 

(P^7^ - M fi)^ (x) = , n = n{xo) = ^a' + h^xl , (2) 
V^ = id,-qA^{x) , [Y,YU = 2r]'^\ r7'^'^ = diag (1,-1, -1,-1), 

where x'^ = r] is conformal time, q is charge of a particle, for example, q = —\e\ for electron. 
The time-independent inner product of the solutions of the equation (H) may be chosen as 

(^,^') = y^(a;)7V(a;)rfx. (3) 

As usual, it is convenient to present ip{x) in the following form 

^P{x) = (P^r + Mn) (j){x) . (4) 

Then the functions (j) have to obey the squared Dirac equation, 
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p2 _ (^^^)2 _ ^^t^up^^ ^ zM9ofi7°j (l>ix) = , (5) 

The external electromagnetic field in our case consists of a constant uniform electric (E) 
and parallel to it magnetic (H) fields, 

Fo3 = i?, = H {5X ~ 6l5l) . (6) 

For such a field we select the following potentials: 

Ao = 0, A3 = Ex\ Ai = = -HX26I 1 = 1,2. (7) 

For 6 = we have a usual flat-space case with the mass m = aM. In this case particle- 
antiparticle classified solutions of the equation and all the GF were found in p2| , p!0| , |TI| . The 



case b ^ is of special interest for us. In the case one can consider spinor field in the 
conformally-flat Universe (with scale factor Q) filled by a constant uniform electromagnetic 
field. Making a standard conformal transformation of the gravitational metric and spinor 
field, we arrive to a theory in fiat space-time with a time-dependent mass (QED-i7). The cor- 
responding field equation is given by (0) (an electromagnetic field should not be transformed 
under the conformal transformation). Note that the such a conformal transformation may 
be used also for interacting theories ||12| , [15| . Thus, Eq. (§) is actually relevant to the quan- 



tum spinor field in the expanding FRW Universe with the external constant electromagnetic 
field. 

Remember [|l^ that the expansion law with a 7^ is necessarily connected with a nonzero 



energy density of the cosmological substratum in the early, radiation-dominated phase of 
the Universe. The minimum value of Q {^min = c^) is caused by the strong interaction. It is 
difficult to find solutions of the equation (|^) with the last spin term. Nevertheless, in many 
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interesting cases with strong background part of this term (for region of small (a^o)^ ~ {(^/^Y) 
may be treated perturbatively starting from the reduced squared Dirac equation: 



where the asymptotic form of the last spin term is used {d^Vl h for large enough xq). If 
the intensity of an electric field is more than the characterizing parameter of a mass-like 
potential, {qE)"^ ^ (bM)"^, or its potential is not intense, a?M/h ^ 1, then both the last 
spin term of equation and of the equation (||) can be disregarded. If {qE)"^ < (6M)^, one 
can not neglect the contribution from the term iMd^VL'-y^ . However, our main interest is to 
study an intense gravitational background when a?M/h ^ 1. Then, as it may be seen from 
explicit solutions of the equation one can calculate any spin contributions to matrix 
elements of operators using solutions of equation (|]) with a relative accuracy of the order 
of a^M/h. And one may use these sets as a basis to construct solutions of the equation 
perturbatively, considering the term iM^d^Vt — 6)7*^ as a perturbation. Moreover, the large 
time ((a;o)^ ^ {(^/^Y ) asymptotic solutions of both equations (|^) and (H) are the same. 

To construct the above mentioned generalized Furry picture for both QED and QED-f2 
one has to find special sets of classified solutions of the equation (Q), namely, two complete 
and orthonormal sets of solution: | ±'^{n}{'^) r ? which describes particles (+) and antiparticles 



and antiparticles (— ) in the final time instant [x^ —>■ +00). According to the general approach 
, which can be easily adapted to QED-fi, such solutions obey the following asymptotic 
conditions 



p2 - {MQ'f 



(8) 



(— ) in the initial time instant 




Ho.pix ) ci^{n}{,x) = , sgn = (, x -00 , 

Ho.p.ix^) ^ip{i}{x) = ^^p{i}{x), sgn = C, a;° ^ +00 , C = ± , 



(9) 
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where (, {n} and (,{1} are complete sets of quantum numbers which characterize solutions 
(ip{n}{x) and ^^/'{;}(a;) respectively, Ho,p, = 7°(Mfi — YVi) is one-particle Dirac Hamil- 
tonian; +£, are particle quasi-energies and |^£:| and \-e\ are antiparticles quasi-energies. 
All the information about the processes of particles scattering and creation by an external 
field (in zeroth order with respect to the radiative corrections) can be extracted from the 
decomposition coefficients, which form the matrices G (cl'' ); 

^V-lx) = +^(a;)G'(+|^) . (10) 

The matrices G (d^ ) obey the following relations, 

G(cr)G(cr)VG(d-)G(d-)' = i, 

G (x) G u^y + G u\-) G (\-y = , (11) 



where I is the identity matrix. All GF in the formalism may be also constructed with the 
sets |±V'{n}(a;)| and |^^/'{„}(a;)|. Below we are going to present such solutions. 
The functions 0(a;) can be written in the following form: 

where Xj_ = {0,x^,x^,0), xj^ = (a;°, 0, 0, x^); {p3,pi,n,C,,r} is a complete set of quantum 
numbers. Among them and pi are momenta of the continuous spectrum, n is an integer 
quantum number, ^ = ±1 and r = ±1 are spin quantum numbers; v^r are some constant 
orthonormal spinors, f^^f^r' = ^rr'- The eq.(||) allows one to subject these spinors to some 
supplementary conditions. 



Ev^r = ^v^r, S = 70(^^7^ - 6M)/p, p = ^{qEf + (6M)2 , 
Rv^r = rv^r, R = sgn (qH) i'j'^'y'^. (13) 
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If if 7^ 0, the function (ppj^nr{x±) has the form 

1/2 



exp |-zpia;^ _ (X) , X = v^|gi/| ( + , 



2"+%2n! 

where Ti.n{x) are Hermite polynomials with integer n = 0,1, . . .. If H = 0, the discrete 
quantum number n has to be replaced by the momentum p2, and the corresponding function 
has the form (ppj^nr{xj_) = (27r)^-'^exp {—i {pix^ + P2x'^)} ■ Let us present the function 
as follows 



(14) 



where (pp^n^rix'^ ,Pz) is a solution of equation 



{Pz - pv f - pX- ipi 



(15) 



pI + pI 



with fj = x° - p-^qEps, pX = pl{hM/pf + uj + a^W^ , 

\qH\{2n+l-r), n = 0,l,... ,H^0 

H = 

One can form two complete sets {±4>p^nir {x'^,Pz)} and {^^pan^r {x'^,Pz)} of the solutions of 
equation ([l5|) using the functions 

1 



\(j)p^n^r [x'^,Pz) = Q/^^-5/2[±(l - i)r], 



r = —{pr] -p, 

'-(Ppsn^r (x°,Pz) = C^^-i/-l+C/2 [±(1 + «)^] , 



(16) 



2 2 

Similar solutions were first presented in p2|. Then, solutions of equation (H) 0(a;) can be 
constructed as follows, 



±4'p3Pin^r{x) ±'Pp3Pin£,r{x, Pz)\p2=0, 
-.(pp3pin^r{x,Pz) = (27r)"^/^e"'P^^^ ±(j)p3n(r{x^,Pz) (l)pinr{x±)v^r 



(17) 
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and in the same form with (±) indices above. 

One can verify that the solutions of the Dirac equation 
with different ^, namely, {V^Y + ±0p3,pi,n,+i,r(a;) and {V^'^^ + MVl) ±(j)p^^p^^n,-i,rix), or 
i'Pfj.l^ + MQ) '^(f)p^^p-^^n,+i,r{x) and (Vfj,'y^ + MQ) ^0p3,pi,n,-i,r(a;) are linearly dependent for 
each sign " +" or " — " . Thus, to construct the complete sets we may use only the following 
sets of solutions: 

(18) 
(19) 

Choosing the coefficients C and C in ([T6|) as follows: C+i = (2p)^^/^ exp(— ttA/S) and 
C+i = ip^)~^^'^ exp (— 7rA/8) , one gets two complete sets { ±V'p3pinr(a;} and {^ipp^p^nrix)} of 
orthonormalized solutions of the equation (Q). These solutions are classified as particles (+) 
and antiparticles (— ) at x° — ±oo according to the asymptotic forms of the corresponding 
quasienergies, = Cpk°| and ''e = (p\x'^\ (see for additional arguments advocating 



such a classification). It matches with classification |^ of similar solutions in QED. 

According to the above discussion the solutions (|18[) and ( [T9| ) of the Dirac equation may 
serve in an intense gravitational background, a?M/h ^ 1, and in the other cases mentioned 
after Eq. (||). Satisfying the Cauchy conditions one can see the solutions (|18]) are vahd with 

< 0, \x^\ ^ a/h and the solutions ([T9| ) are valid with > 0, ^ a/b for any intensity 
of the background. 

To find the matrices G (cl^') defined by (|10D, it is convenient to use an asymptotic form 
of the solutions. Using (H) and (|13D we get 

G = 6i,, g , / = (p3,Pi,n,r), I' = {p-,' , p,' , n' , r') , (20) 

where 
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gicf) = cCn,+l,r- ^ ^0 {ido - PV) ^' (f)p,,n,+l,r {x'^,Pz) ■ (21) 

III. GREEN FUNCTIONS 

Let us start with out-in GF which is the causal propagator 

S'^{x, x') = c~'^i < 0, out\T'ip{x)'ijj{x')\0,in >, =< 0^ out\0,in > . (22) 

Here ipi^) is quantum spinor field satisfying the Dirac equation (H), \0,in > and \0,out > 
are the initial and the final vacuum, and Cy is the vacuum-to-vacuum transition amplitude. 
The propagator S'^{x, x') obeys the equation 

(P^7^ - Mfi) S\x, x') = -6^^\x - x') . (23) 

Another important singular function is the commutation function S{x,x') = 
i tp{x),ilj{x') ^. It obeys the homogeneous Dirac equation (^ and the initial condition 
S{x,x')\^^^^, = i7°5(x — x'). Besides, the following GF are studied 



S^n{x,x') = i< 0,in\Tip{x)ip{x')\0,in >, S1^{x,x') = i< 0, m|'^(x)'?/'(x')T|0, in >, 
S"" (x, x) = i < 0, in\ip{x)'il){x')\Q, in >, S^{x, x') = i < 0, in\ip{x')ip{x)\0, in >, 
S^^t{x,x') = i < 0, out\T'ilj{x)ip{x')\0, out > . (24) 

Here the T-product acts on both sides: it orders the field operators to the right of its and 
antiorders them to the left. The functions Sf^ and S^^^ obey the equation (|^), 5*^ satisfy 
the equation (j|), and S'[^ obeys the equation {V^^^ — Mil) S^^{x, x') = 6^'^\x — x') . 
One can express the GF via the solutions (p^Sf) and (^ 
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{x, X' 
S (x, x 



^e{xo- 4) s- (x, x') - ^ (4 - xo) (x, x') , 

= 8- {x,x') + S+ {x,x'), 

^e{xo- 4) 5- (x, x') - ^ - X,) SI {x, x') , 
= ^ K - a^o) -^i; (2^, a;') -e{xQ- x'q) 5+ (x, a;') , 

^e{xQ- x'q) s-^i {x, x') -e{x'Q- xq) 5+ 1 {x, x') , 



(25) 



where 



-oo nr 

/~^°° r ll * ^ 

dpsdpiY^ -i'pspmrix) g{-\~y ~i>p3pinr{x') , 
-oo 

/+00 ^ 
dp3dpi X, iV'pspmr (a;)± ^pspmr {x') , 
-oo jjj. 

/+00 _ 
dp3dpiJ2 ^fPp3Pinr{x) ^Ippspinr {x') . (26) 
-oo nr 

J2nr means the summation over all discrete quantum numbers n, r (and the integration over 
the continuous p2 ii H = 0). Using the relations between GF and between the matrices 
^ (cl^ ) present the functions 5"^, and S'^^ as follows 

±S^{x, x') = S%x, x') ± e{^{xo - x'q))S{x, x'), 

^S^nix, x ) = SIX^', x/) ± 9{t{xo - Xq))S{x, x'), 
±S^,{x, x') = SUx: x') ± e{T{xo - x',))S{x, x'), 

^in{X} X ) = S {x, x ) — S {x, x ), S^^i{x, x ) = S {x, x ) — 5'^(x, x ), 



dpsdpiY^ -i/jpspinrix) g{+\~)g{- 

-oo „^ 



-oo 



S^{x,x')=i dpsdpij^ ^'^P3Pinr{x) g{+\'^) ^g{+ 



+'^P3Pinr(-^ )■> 
) ^pspmrix ) 



(27) 



12 



Let us consider the functions and S"'''^. The coefficients ( pi]) do not depend on pz, 
thus, one can present the functions and S""'^ in the following convenient form 



x') 



hoo /■+00 

(i; 



dps ^-ipsy^^ 



27r 



) i/M -^A* "^A* ' 



(28) 



where 

Qa 

and 



ZiPz)9 (^ + \ ) +'^P3Pinr {j]'}X±} Z ,Pz^ , 

nr 

dpzdpiJ2 -'ipp3Pinrifj,X±,Z,Pz) [^(-1")" J ~^p3pinr , 
-oo „^ 

dPzdpiY^ -^P3P^nr{V,X^,Z,Pz) [g{+\')9{-\'y \ +tp p^p^^^iv' , x'^, z' , Pz) , 
i\ dpzdpiY, ^^P3Pinr{V,X±,Z,Pz) gi+l^)' 9 ' i^p,p,nriv' , x'^, z' , Pz) , (29) 



± 



^pspmr (^?, 2,^2) = (7°zi9o + 7^ (Ps - qEx°) + 7_l (z9 - qA) + MQj^(f)p 

(x_L)t'+l,r, ±4>p3nr{'n,Z,Pz) = +i,r(a; ,^2), (30) 



27r 

and in the same form with (±) indices above. Within the same approximation one can 
rewrite the functions SS'"'^ as follows: 



S^'"'" = {I'tdo + 7' (P3 - qEx^) + 7x {^^ - qA) + M^) A]^'^ 



where the functions Ag "'^ obey the equation 



a,p 



(31) 



d\ ( .d X . plibM^] 



q 



A 



0. 



They may be easily expressed via the solutions ( |30D according to the Eqs. (^). The functions 
Aq'"'^ are just the GF of the squared Dirac equation in electromagnetic background p!0| , p 
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where f] is the time, z is the coordinate along the electric field, the mass tuq = p\{hMY / p^, 
the potential of the electromagnetic field is = fjp/q, and the spin term qE'j^'-f^ is changed 
to pS. In what follows we are going to use the following representations ||11|| : 



^Q= fQds , Aq = sgn(?/o) / fqds 
A^= / fQds + 9iz'-z) I fgds, 

JTa •^F3 + F2-Fa 



A^ 



fQds + 9{z-z') / fqds 

Fa "'F3 + F2-Fa 



(32) 



where 6{0) = 1/2, the contours of the integration are indicated on the Fig.l, and 





Im s 

Tc Res 
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1 ► » 

r2 TT 




^2p 


-* 1 





FIG. 1. Contours of integration Fi, Fa, Fc, Fa 



e^''^exp(--gF>^'^.)/< 



(0) 



(33) 



= exp IJ A^dx^^ 4 (z - z') h, 
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f± = (4vr)- 



qH 



sm{qHs 
P 



exp |-^y_Lg-^coth(gFs)?/_L| 



■ P 



P 



exp < —i-{ri + T]')z ~ impS + i— z — {rj — r]') coth(ps) 



sinh(ps) """" I "2' ' ■ "■ "^^ ' "4 

One can calculate in ( |28D all Gaussian integrals over ps and z. Thus, one gets 

S^-^ {x, x') = {YV^ + MVL) A(-) (a;, x') , 



A'^(x, x')= / f{x,x',s)ds, A(x, x') = sgn(?/o) / f{x,x',s)ds 



A"(x,x') = -A«(x,x') - A(2)(x,x'), 
AP(x,x') = -A«(x,x') + A(2)(x,x'), 

/ "'rg+ra+r^ 



A(2)(a;,s') 



fr{x,x',s)ds, 



/rs+ra-r, 

where fr{x,x',s) = 7r^^/^/2 7 (1/2, a) /(x, x', s), 7 (1/2, a) is the incomplete 
function, and a = e'^'^l'^ [^s(^MfujY^ [(xq + x'o) si^Mf + qEy^f . Here 

/(x,x',s) = exp (pSs -^a'^^Fj-^s) /^^^ (a;,x',s) , 



/(°)(a^, a;', = exp A^rfx^ j /y /. 

. P /-^^ ^ , ' ^ 3 P 



exp<{2^(xo + x'o)l/''-«7(a;o-x'o)^coth(ps) 
smh ( ps ) (-J ^/"^ 12 4 



-i(aM)^s + i-^yl coth(ps) - {hMfs (xq + x'o)^ + 2gE?/^ (xq + x'o 
4a; AlOJ l 

where u; = s coth(ps)(6M)^/p + {qEY / p^. One can see that 



d 



/(x, x', s) = (^M2^]2 _ p2 ^ I^^^^F^. - ^6M7° ) /(x, x', s) 
lim /(x, x', s) = i^'-^Vx — x') . 
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Thus, f{x,x',s) is the Fock-Schwinger function PJ2^ of the QED-f2 theory, and s is the 
Fock-Schwinger proper-time in this case. The contour — r2 — Fi in representation (^) for 
A transformed into F (see Fig.2) after the integration over and z. The resuh is consistent 



with the general expression for the commutation function obtained in The function 

f^'^^ (x, x', s) which has appeared in (^) coincides with the Fock-Schwinger function of the 
scalar case. Due to that we can widely use results presented in 

Im s 



O 



Re s 



FIG. 2. Contour of integration F 

If 6 7^ 0, then the function f{x,x',s) has three singular points on the complex region 
between contours F^ — Fi and F^ — F3, which are distributed at the imaginary axis: pso = 0, 
psi = —m and ps2 = —ic2- The latter point is connected with zero value of the function uj. 
We get an equation for C2 from the condition u; = 0, C2tan(c2 — 7r/2) — {qE/{hM)Y = , 
where 7r/2 < C2 < vr. The position of this point depends on the ratio qE/{bM), e.g. C2 — > vr 
as bM/ {qE) — > and C2 ^ 7r/2 as qE/{bM) 0. Notice, in the case ii^ = it is convenient 
to put C2 = vr/2 + since the contour F2 is also passed above the singular point S2. If 6 = 0, 
then UJ = 1, and the function /(x, x', s) has only two singular points sq and Si on the above 



mentioned complex region. In this degenerate case we have known representation |11 . 
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Similarly to Ref. one may demonstrate that the function A^^^ from ( ^5]) can be pre- 
sented via proper-time integrals with the Fock-Schwinger kernel f{x, x', s). Then, using (p7|), 
one can verify that the representations ( p5|) for A*^^^ and A^^^ obey the equation (|^). Thus, 
all the A-functions considered here, excluding those marked by the index "c" , are solutions 
of the equation @. The important difference between the functions A'^, A^^^ and A, A*^^) is 
that the first ones are symmetric under simultaneous change of sign in xq, x'q, x^, x'^ and 
the second ones change sign in this case. 

Using the kernel (^) one can express out-in effective action in the QED-fi theory, 

r roo 

^out-in =2/^^/0 "^"^ "^^^^ ^' ^^^^ 

IV. VACUUM INSTABILITY, MEAN VALUES OF CURRENT, AND ENERGY 

MOMENTUM TENSOR. DISCUSSION 

Using the exact solutions and the GF constructed above, one may calculate different 
physical effects having both local and global character. The proper-time representation of 
the GF gives us a possibility to study all mean values of current and energy momentum 
tensor in the same manner. 

All the information about the processes of particles creation, annihilation, and scatter- 
ing in an external field (without radiative corrections) can be extracted from the matrices 
G{cf) (0). These matrices define a canonical transformation between in and out creation 
and annihilation operators in the generalized Furry representation 

a\out) = a\m)G (+|+) + b{in)G (_|+) , b{out) = a\in)G (+r) + b{in)G (_r) . (40) 
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Here a] (in), h\{in), ai{in), bi{in) are creation and annihilation operators of in-particles and 
antiparticles respectively and al{out),bl{out), ai{out), hi{out) are ones of out-particles and 
antiparticles, / are possible quantum numbers (in our case it / = pi,p3,n,r) . For example, 
the mean differential number of particles created (which are also equal to the number of 
pairs created) by the external field from the in- vacuum |0,m > is 



Ni =< 0, in\al{out)ai{out)\0, in >- 



9[- 



(41) 



(for a review of gravitational particles creation, see |2^J2^). The standard space coordinate 
volume regularization was used to get the latter formula, so that 6{pj — p'j) ^pj,p'-- The 
probability for a vacuum to remain a vacuum is 



P, = |c,p = exp|^ln(l-iV,)| 



(42) 



Similar to the electric field case 21 we get 



N, 



-it\ 



if > 1, and ^pT > A, and p^T > \qEp 



'31 , 



(43) 



where A is defined in (0). The latter conditions take place with large enough time for 
action of the electric-like field, T = — a;?„ . The result (^) coincides with one obtained 



m 



I9|, and for 6 = it coincides with one obtained in [pT|. Evidently the creation process 



is a coherent effect of both the electromagnetic and gravitational fields. If the condition 
p\{hMY / « 1 takes place (the gravitational field is in a sense weaker than the electric 
one), the ps dependence on Ni is similar to the case 6 = 0. Thus , one can estimate that 
/ dp^ = {qE)~^p'^T |^T|. Then the particle creation per unit of time may be calculated 
similar to ||18|. In strong enough gravitational background the time dependence of the effect 
is nonlinear and needs to be studied specially. 
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To get the total number N of particles created one has to sum over the quantum numbers 
/. The sum over the momenta can be easily transformed into an integral. Thus, if 6 = one 
gets result presented in [^. If 6 7^ 0, the total number of pairs created per space coordinate 
volume has the form 

,,=. = ^ = ^!m^expj-.(^|, (44) 

/ cix 47r2 6M \ P J 

where (3{n) = qH coth{mTqH/ p). The observable number density of the created pairs in 
the asymptotic region xq = Xg"* ^ 00 is given by the expression n'^^ = n'^^ /Q^{xo) ■ If the 
electromagnetic field is absent and a = , these results coincide with ones in |^D[. In case 
6 — >■ the expression (^) is growing infinitely. In this case the particles are created in main 
by the electric field, whereas the parameter b plays a role of "cut-off" factor, which eliminates 
creation of particles with extremely high momenta along the electric field. It is seen from 
the expression (|43|) . Thus, the limit 6^0 corresponds to the case of the electric field which 
acts for infinite time. Then the number of particles created is proportional to the time of 
the field action. As was already remarked above, in this case / dp^ = {qE)^^p'^T. Then, 
the parameter b may be understood as (y^)(TM)~^. The vacuum-to-vacuum transition 
probability can be calculated, using formula p^). Thus, we get an analog of the well-known 
Schwinger formula in the case under consideration, 

P„ = exp{-Kn"/rfx}, ^=£^-^(^±|^exp{-„.<i^^} . (45) 

Now we are going to discuss vacuum matrix elements of current and of metric energy- 
momentum tensor (EMT) of spinor field. Making the conformal transformation the current 
in FRW universe may be presented in the following form 

J, = n-'{x')j^, J, = I [^(x), ^(x)] , (46) 
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where is QED-fi current of the spinor field ip{x). The latter obeys the Dirac equation 
(0). The EMT of spinor field in FRW universe may be presented as 



T, 



rpcan 

flU 



ncan ^ 

■ VjJ. I 1 1 



(47) 



where T^^i, is EMT in QED-f2. In the case of unstable vacuum one needs to calculate three 
types of matrix elements |^ , depending on the problem in question: 

< 0, out\ji,\0, in > c'^ , < T^u >"=< 0, out\T^^\0, in > , (48) 

< 0,m|j^|0,m > , < T,,, 0,m|r^,|0,m > , (49) 

< >°«*=< 0, out\j^\0, out> , < T^^ >""*=< 0, out\Tf,^\0, out > , (50) 

Using GF, which were found before, one can present these matrix elements as follows: 



< tqti {^^S'^ix, x)} = iq tr {7^ {^^V. + Mfi) A'{x, x')}^^^, 

< T,, z/4 tr { (7^ {v. + V'l) + 7. {V, + V'l)) S\x, x')}|_^, 
= iti{B^,/\%x,x')}l^^, , 

< 3, J, >' + < J, >^'^ + < J, , 



< J, , 



^ rri ~^OUt_^ rp -^C , ^ rp (1) _ ^ rp (2) 

< >(i'2)= tq tr {7^ (7'^P, + Mfi) A^'^^\x, x')}l^^, , 

< T^. >(i'2)= ^ tr {i?^,A(i'2)(x,x')}[^^, , 

B,, = 1/4 {7^ {v. + V'l) + 7. {V, + V'l) } (7'^P. + Mfi) 



(51) 

(52) 
(53) 
(54) 
(55) 
(56) 
(57) 
(58) 
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where P'* = —i-Q^rp: — Q A^{x'), the GF are given by Eqs. ( (|35|) , and the relation A'^(x, x) = 
(1/2) [A'{x,x) - A+{x,x)] is used. 



The contour r2 does not pass through the singular points of the function f{x,x',s), thus 
the contributions from A^^ to EMT and to the current are finite. Displacing the contour r2 
to the real axis we may present the difference A'^ — A^^ as 



where A'^ is related to the anticausal GF, S'^{x,x') = {jV + MQ)A'^{x,x'), and the integral 
A'" can be expressed by means of the A-function related to the commutation function, 
S{x,x') = {jV + MQ)A{x,x') A^{x,x') = sgn(xo — x'q)A{x,x'). Taking into account 
the above definition the space-like part of the limit x — x' ^ in the EMT expressions 
have to be treated either as the limit xq — Xq = +0 or as xq — x'q = —0, respectively. Then, 
in according to the initial condition for the commutation function we get A{x,x')\^^^^, = 
0, doA{x, x')\^^^^,^ = 5(x — x'). Thus, the contributions from A into < T^^ and < >*"■ 

are zero. (All contributions from A, which might appear as a result of a change of the limit 
definition, are background independent. They may be eliminated by a renormalization.) It 
follows from (^) that 



To get convenient forms of < and < T^^, >*" we may rewrite A*^^^ as follows: 




(59) 



A^ix, x') - A^^x, x') = A%x, x') + A^(x, x') 




(60) 



f{x,x' 



s*) = i'f{x',x,syj'. 



(61) 
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Changing s — > — s one can represent integral from (^) in the form 

poo 

A%x,x') = - -f'^f{x',x,sy-f°ds, (62) 
Jo 

to get 

< >^"= Re < >' + < >(2) + < >(3) , (63) 

< T^, Re < T^, >^ + < T^, >(2) + < T^, >(3) , (64) 

where the terms < >^^^ and < T^^, >'-'^'' represent the contributions from A*^^-*: 

<j^>(3)=,gtr{7^(7''P. + Mn)A(3)(x,a;')}[^^, , (65) 
<T^, >(3)=ztr{5^,A(3)(x,x')}L^^, . (66) 

Using the Eq. (|6T|) one can verify that each of the terms < >*^^\ < >^^\ < T^y >(^) 
and < Tfj_y >^^) is real, as it follows from their initial definitions. 

One can get similar expressions for the out-out-matrix elements of the current and EMT: 

< 3, = Re < 3^. - < J, >^'^ + < >^'^ , (67) 

< T^, Re < T^, >^ -< T^, + < T^, . (68) 

All nondiagonal matrix elements of EMT are zero: 

< >^=< T^, >(2)=< T^, >(3)= 0, fi^u. (69) 

Using the formulas 

exp (pSs) = cosh (ps) -|- S sinh (ps) , 

1 

-( 
2 



exp ( -ij^a^^'F^^s ] = cos (qHs) + -f^-f^ sin (giJs) 
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one may calculate the traces for the current and EMT: 

t(s) = trF = 4 cosh (ps) cos (qHs) , Y = jexp (^pEs - i^a'^'F^^s^ | , 
tr{Sy} = tanh(ps)r(s) , tr{7V^} = ^ tanh(ps)r(s) , 
tr{7°y} = -^tanh(ps)r(s) , tr{7^7^r} = -tSin{qHs)T{s) , 

and other traces are zero. Then, the current matrix elements can be found in the forms 

<Jm>'=/ a^{s)r{s)f^''\x,x',s)ds\ (70) 

JVc \x=x' 



-- / ai,{s)T{s)fj.°\x,x',s)ds 

Jr3+T2-Ta 



(71) 
(72) 



qE 

V3 + Vo-^teinh{ps) 



where the only x^ components (along the electric field) of the currents differ from zero and 
vanish in the absence of the electric field, as it presented to be. The nonzero components of 



< T^^ >^ are 



< 



T^.>''= I t^{s)r{s)f'\x,x,s)ds, iip^v 

qH 



(73) 



tJs 



sinh(2ps)' sm{2qHs) ' 

2 ^qE 



t3(s) = Acoth(ps) - ^^^^ 



— (a;o + x'q){1 - uj'-^) 



200 ' 2p 

and nonzero components of < 7)^^ >(^) and < T^i, >^^^ can be written as follows: 
<T^,>i^)^-\[ l^{s)r{s)f^''\x,x,s)ds, iifi^u, 
< T^u / l,{s)T{s)f^'\x,x',s)ds , , iill^U, 



(74) 
(75) 



li{s) = ti(s), ^2(5) = t2{s), his) = i {Vsf - ■^^tanh(ps) , 



/o(s) = h{s) + his) + his) + i 



2p 

2 {uj-^h^M'^sxof - iuj-^b'^M^s + M^Q' 



- tanh(ps) . 
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Here, lo{s) is obtained from to{s), using Eq. (^) and taking into account that the A'^^) 
function obeys the equation (H). 

The components < and < T^j, are expressed by the causal GF in the well- 

known Schwinger's form. They are sources of information about local features of the theory. 
For the first time it is obtained exactly with respect to the given external background. 
The components < >(2>3) ^ j^^^ >(2,3) ^^.q Qj^jy related to global features of the 
theory and characterize the vacuum instability. Such expressions cannot be calculated in 
the frame of the perturbation theory with respect to the external background or in the frame 
of the WKB method. Then, such terms has never been studied. If the parameter b 0, 
the only expression (^) for < T^^ >^ has to be regularized and then renormalized. The 
expression (^) for < is finite after the regularization lifting. The terms < >(2.3) 

and < T^^ >(2:3) g^j^g g^jgQ gnite. Thus, as one can see from Eqs. (^) and (^), the divergence 
in < T^i, >*"■ and < T^^, >°"* is the same as in < T^j,^ . That is consistent with the fact 
that the ultraviolet divergences have a local nature and result (as in the theory without an 
external field) from the leading local terms at s — +0. Nevertheless, as 6 — > ( flat space 
limit) the terms < >(2.3) g^^^j ^ j-^^ >(2,3) y^q^q electric field dependent divergences defined 
by the dimensionless parameter f3 = bM/qE. The nature of that is similar to the divergence 
of n"' for 6^0. In the case 6 = the density of excited vacuum states is proportional 
to the time T. After a special regularization with respect to time T all of the terms 
become finite. 

Now we are going to discuss current, energy density and pressure of the created particles. 
Let us introduce normalized values of current and EMT, which may be easily connected to 
observable values in some appropriate asymptotic region xq = Xg'*, 

f; = ~f;/^'M , = f-/n\x,) , (76) 
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where according to the definitions p9|) and (0) the corresponding densities of particles 
created per space-coordinates volume are 



^i^u- J^^^ 5 Xq-Xq . {(^) 



Then, using representations (|53D, ( |51D and (§3), one gets from ( [77D and ( [75[ ) 



j7 = 2<j;,>(2), f- = 2<T^, >(2), xo = xS^ (79) 

The component Tq^ ° = f2~^TQQ is the energy density and —T^ ^ = f2~^r^J^ for /i = z/ = 
1,2,3 are components of pressure which are measured by the cosmic observer relative to the 
measured volume. In contrast to that, measured energy and pressure taken per coordinate 
volume are given by f^-" ° = ^-^f^^, -f^'^ " = Vt-^f^l, /i = z/ = 1, 2, 3. 

We are going to analyze contributions to the quantities (^) at xq = Xq*. The leading 
asymptotics in < >'^^) and < T^y >(^) at a;o >> ^ji^M) are determinated by the 
expressions 

3, >^'^= Ms,) tr{A(2)(a;,x')}[^^, , T,, >(2)= ^s,) tr {a(^)(x, . (80) 

An asymptotic expression for trA'-^^ with x ^ x' can be found using the method described 
in the App. A of Ref. [|l| for the scalar case. Then, one gets 

where n'^^ is defined in (^4]). This expression is also valid when xq < 0, since the function 
A*^^) is odd in xq as x = x' (see Sec. Ill ). Thus, one can see when xq >> y/p/{bM) the 
leading terms in < >^'^^ and < T^j, >^^'' are 
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< Too >^^^= [p'^xl + a^M^ + 2qHsmh-^ i^nqH/p) 

< Tn >(^)=< T22 >^^^= qH sinh-^ i^nqH/p) 

< T33 >(^) 



pxo ' 



(gExo)' + ^ 



pxo 

3 

lb , . 

— ■ 81 

pXo 



2Ti{hMf 

Comparing our results with the scalar theory case |T[|, when xo c)o we may see the 
difference only in the factor n^*" which enters in the leading terms of < >*^^-', < Too >''^'* 
and < T33 >(2) . 

Doubling the expressions (0) according to the Eqs. (|79D , one gets the mean densities 
for current and EMT of particles created. It turns out that these quantities are proportional 
to the density of total number of particles and antiparticles created (2n'^'') for the infinite 
time and do not change their structure with increasing of xq. The latter means that one can 
consider all the expressions obtained at any fixed Xq as asymptotic forms if Xq ^ ^jhM. 
In a strong background a^M"^/ p < 1 and, therefore, such a time has to obey the condition 
Xo >> a/h. Thus, in the strong background our asymptotic conformal time Xq, which is 
large enough in quantum sense explained, corresponds to the large cosmological time t. 

Note, that one can neglect the second term in the brackets of the expression ( |8lD for 
< T33 >(2) at hM/{qE) < 1. Also one can neglect both the term a?M'^ in the brackets of 
the expression (^Tj) for < Too >*^^'* in case of strong background a^M^/p < 1 and third term 
in the same expression if the magnetic field is not strong enough, qH/p < 1. The current 
density j^'^ = 2 < >*^^^ does not depend on the asymptotic time. As 6 ^ this expression 
coincides with the one for flat space, < >^'^= —2qn'^^6^. The pressure component along 
the electric field direction T3"' = 2 < T33 >'^^) is growing with time upon the action of the 
field. However, if qE/{bM) << 1 then the asymptotic condition for xo is consistent with 
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the fact that the term [qEx^Y the expressions < T33 >(^) from (|8lD will not be dominant 
before large enough time instant xq. In this case one can neglect the contribution which 
depends on the electric field, if the field is switched off at a fixed time. 

The components of the pressure in the directions, which are perpendicular to the electric 
field, = 2 < Tn >(2) and T!^!^ = 2 < T22 >^'^\ decrease when the magnetic field increases. 
If an electromagnetic field is absent, the pressure is isotropic according to the symmetry 
of the space-time and to the corresponding symmetry of the vacuum. For 00 the 

remaining terms of the measured energy and pressure of the created particles taken per 
coordinate volume are 

rpcr r)~cr rpcr S 

{qEf 

n~cr 

1, - ^Zn , 1, - Zn , 

and, if an electric field is absent the only remaining term is Tq^ ° = 2Mff'^ . The last repre- 



sents the total rest mass per coordinate volume and coincides with the result of Ref. [|T8 . 

Another problem is to study a back-reaction of particles created on the electromagnetic 
field and metrics, or to be more correct, it is better to say about back-reaction effects pro- 
duced by both of particles created from a vacuum and polarization of an unstable vacuum. 
To this end one needs to use the expressions < and < T^^, >*" for all times Xq. The 

explicit proper-time expressions found above give us promising tool for accurate analysis of 
the total back-reaction related to the vacuum instability. In this approach we do not need to 
select phenomenologically parts come from real particles and from a vacuum polarization, 
that has been usually done in literature (see, for example,, review in [^). Of course, to find 
self-consistent solution taking into account the back-reaction one needs numerical estima- 
tions (compare with pure electromagnetic case, [^). Keeping in mind such an application, 
one needs to get preliminary information about the behavior of the expressions (|63|) and 
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(|6^) in time and to select the leading components. To this end let us estimate the above 
mentioned expressions at characteristic large time, » p/{bMY, and at characteristic 
small time, Xg << p/ibM)"^. 

Neglecting divergent terms in < T^j, according to the standard renormalization pro- 
cedure, one can reduce Eq. ( p^ ) to the following finite form 

< T^, Re < T^, >),^ + < T^, >(2) + < T^, >(3) . (82) 

Here < T^^ >}in is a finite part of < T^u >'^. An estimation of < T^,, >}in can be made 
only after renormalization. Since we are here interested in to reveal global features of the 
theory, details of the renormalization problem together with others, related to the renormal- 
ization, will be considered in the next paper. To select contributions related to the vacuum 
instability, we note that the functions (^) and A*^^-* ( ^91) with x = x' are even in xq. 
Thus, the functions < T^u >^ and < T^u >^^^ are also even and do not vanish when xq 0, 
whereas the functions < >^ and < >^^^ are odd and vanish in the limit. Moreover, 
< >^=< >(^)= for all Xq with b = 0. The proper-time integral A*^^) ( |55D is odd in 
Xq with X = x' and vanishes when xq 0. Thus, the expression < T^,^ >*^^'' is also odd in 
Xq and vanishes in this limit. The term < >^'^^ is not zero ii E 0. It is even in xq and 
differs from zero when xq — > 0. 

The asymptotic behavior of < >^^^ and < >^^^ is defined by the asymptotic 
expression for A*^^). Using the method presented in Appendix A of Ref. [|l| one can verify 
that asymptotically 

A(3)(x,x') = sign(xo)A(2)(x,x') , (83) 

and therefore. 
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< sign(xo) < J;. , < T^. >(3)= sign(xo) < T^. >('K (84) 

The expression ([7^) does not need to be renormahzed, thus, one can easily verify that the 
relation < >'^~ Xq ^ — > holds asymptotically. Then, the asymptotics of < is 
determined by the asymptotic behavior of < >^'^^ and < >'''^'* . If Xq with the large 
modulus is negative, < >™= Re < and < T^i, >Yin= < ^^i^ ^/m- Thus, in 

this regime the vacuum instability in the time-dependent background does not affect the 
expressions. If /5 = bM/{qE) « 1, the expression Re < T^j, >jj„ remains finite with 
6^0, while the terms < T^,^ >^'^^ and < T^^ >^^^ diverge in the limit since the global effect: 
the total density of excited vacuum states in an electric field grows infinitely and its mean 
momentum along the electric field increases. In this case the asymptotics of < T^i, >Yin 
determined by asymptotic behavior of < T^,y >*^^^ and < T^,y >^^\ i.e., by the asymptotic 
global properties of the theory. 

Let us estimate < >*" and < T^^ >j"„ for small xq. The terms < >^ and < >^^^ 
are odd in xq. They vanish at Xq — ^ 0. That is why the leading term in ( |53D is < >^'^^ . 
The leading contribution from the latter can be represent as 

<j,>^'^=a,{s2)ti{A^'\x,x')} , . 

1. J x=x' 

Then, using the expressions (|A1|) from the Appendix A, one gets 



l + (c2{hMf/{pqE))' 



€ (85) 



The expression (^) differs essentially from the asymptotic one (0). That remains true in 
the quasi-flat metrics ( f3 = bM/{qE) << 1), when the vacuum instability results only from 
the electric field. In the case of the small time limit the expression < >(2)= qff^'5^^ differs 



by a sign from the asymptotics (|8l|). The factor n'^^ is known from the asymptotics. That 

29 



is natural since value of the time instant Xq is large enough (xq + T/2 >> (gi?)"^/^ and 
Xq + T/2 >> (aMy(qE)~^^'^ according to (^31) ) with respect to the initial time instant 
(xq" = —T/2). Then the field action time is large enough to obey the stabilization condition 
and the density of exited states has already the asymptotic form. However, such xq is 
still less than the asymptotic time for any f3. Hence, at the time instant Xq the vacuum 
differs essentially from \0,out > vacuum. Technically that means that at the small time, 
Xq << y/qE /{bM), one cannot obtain the normal form of operators by means of the A°"*- 
function and, consequently, to use only the term < >*^^^ for calculation of the mean current 
of quasiparticles, which are created up to that time. The small time expressions one can use 
only to calculate the back-reaction. Thus, in the external background under consideration a 
small time back-reaction to the mean electromagnetic field is not a screening (as it may be 
expected by analogy with a back-reaction of real particles created) but an induction. Since 
the small time considered is related to early stages of the Universe, such an observation may 
be important. 

Since the function < T^j, >*^^) is odd in xq and vanishes as Xq 0, the leading contribu- 
tion to dH) is determined by Re < >jj„ and < >(^) . However, ii (3 = bM/{qE) « 1, 
the situation is more simple in the domain where >> qEx^, and the leading contribu- 
tions come from < T^i, >^^^ only. That is related to the same global effect, mentioned above 
in the asymptotic case, since the expression Re < T^,y >'fin remains finite with 6 — > 0, while 
the term < T^j, >*^'^^ diverges in the limit. The total density of excited vacuum states grows 
infinitely and the mean momentum of the states along the electric field increases, as well. 
Taking into account that the leading contributions in < T^^j >^^^ result from the integral 
over a neighborhood of the singular point Si, one gets 
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< Too >^')=< T33 >(^)= , 

up 

< Tn >(^)=< T22 >(^^= qH/^smh^^{2TiH/E)n^'l3\nl3-^ . (86) 

We see that leading vacuum polarization effects are in main defined by the total density 
of excited vacuum states, which depend on the complete history of a state. Thus, it is a 
nonlocal contribution. 

Considering the small time limit of < T^^ >'^^) we get 

<T,, >(^)=/^(s2)tr{A(2)(x,a;')}L^^, , 
where trA*-^^ was found in (|A1|) . If /3 << 1, the result has a simple form: 



< Too >^'^=< T33 qE 



< Tn >(2)=< T22 >(^)= -2qHsiYi\i-\2TiH/E)xoTi(i 



2~cr 



One can see that expressions for < T^y >^'^^ at the asymptotic time and at the small time 
are quite different. The only components < Too >''^'*=< T33 >*^^) for /3 << 1 coincide. 

The technics developed in the article allows one take into account global effects of both a 
real particles creation and a polarization on the same footing. The explicit form and the limit 
estimations for the mean vacuum values of EMT may be used for the calculations of the back 
reaction from the unstable vacuum to the gravitational background (in particular computing 



simulation similar to ||3^, see also and references therein). We see a behaviour of such 
components in time are quite different, and one needs to take into account characteristic 
polarization effect. The proper-time representations of the GF may be the necessary step 
in the study of chiral symmetry breaking in QED and the four-fermion models under the 
action of gravitational and electromagnetic fields . Such a study may have an immediate 
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important application to EU, for example, through the construction of inflationary Universe 
where role of inflaton is played by the condensate < ^^\E' >. One can also analyse symmetry 
breaking phenomenon under the combined action of these fields in the Standard Model 



(using also its gauged NJL form [0), or GUT in the same way as it has been done in 
curved spacetime (without electromagnetic field) |T2[. Our methods maybe applied to the 
study of SUSY NJL model in curved spacetime introduced in The chiral symmetry 

breaking under the action of weak gravitational field and constant electromagnetic field has 



been studied in |]3J]. Having our results we may study the chiral symmetry breaking for 
such model in FRW Universe with constant electromagnetic field. 
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APPENDIX A: SMALL TIME EXPANSION OF A(2)-FUNCTI0N 



To calculate small time expansion of tr given by Eq. (pSj) , in the case 



xi « 



p/{hMY and x —>■ x', one can use the method presented in App. B of Ref. [jll. The final 
result is 



tr{A(2)(x,x')} = e'^^{ 
(6M)4 



+ t 



UqEp 



^K{2) + {l/2){xo + x',){y:,fqEK{Q) 



(Al) 
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X exp ( -y^qF coi{TX qF/p)y^ j , 



4+1 (p/(6M))^ 



, l^mvp - c2-i(/ - 1/2) - f 



-C2 j + 2g/7/p sinh-i(2c2gif/p) + ^cl^^ 
n(2) = -l/2r(s2)n(2) , 

(2) 





2 




'4 + {m 


4" 


cl + 


\bM ) 




4 



(2) 



1 












VfcAf ) 


{hMf sm{c2qH/p) 



If bM/qE << 1, the coefficients K(l) and n^"^^ from have more simple form. Thus, in 
the case of an intensive electric field {a^M"^ /{qE) < 1, \H / E\ < 1) we get 



n 



(2) _ ^cr 



l~cr 



(A2) 



Since the final formula (B4) in App. B of Ref. [jT|, which describes small time expansion of 
A^^^ in scalar case, was written with some misprints, we would like to represent its corrected 
form here: 

AfJ{x,x') = e''^''{[z{xo + x'o)c2{bMf/p-qEy'] {-l/2)nfj /{qE)^o 

mr 



+ I 



Ksc{l) 



[Xo + Xq 



l2qEp'^ 
cT {p/{bM) f 



Ksc{2) + (1/2) (xo + x'^){y;fqEK,,{Q) 
, [a^M^/p - c^\l - 1/2) - c-^^ ( 

ci^m^m^ ^^^^ 



+ {qH/p) coth{c2qH/p) + 20^^ 



r2 

'-2 



^ \bM J ^ \bM J 



n 



(2) 
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